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The Fermi-hAT has revealed that rotation powered millisecond pulsars (MSPs) are a major contributor to the 
Galactic 7-ray source population. We discuss the 7-ray emission process within the context of the outer gap 
accelerator model, and use a Monte-Calro method to simulate the Galactic population of the 7-ray emitting 
MSPs. We find that the outer gap accelerator controlled by the magnetic pair-creation process is preferable in 
explaining the possible correlation between the 7-ray luminosity and the spin down power. Our Monte-Calro 
simulation implies that most of the 7-ray emitting MSPs are radio quiet in the present sensitivity of the radio 
survey, indicating that most of the 7-ray MSPs have been unidentified. We argue that the Galactic Fermi 
unidentified sources located at high latitudes should be dominated by MSPs, whereas the sources in the galactic 
plane are dominated by radio-quiet canonical pulsars. 



1. Introduction 

The Fermi Large Area Telescope (Fermi-JjAT) has 
discovered over 80 7-ray pulsars, and revealed that 
the 7-ray pulsars are a major class of Galactic 7-ray 
sources. Of these, Fermi-LAT first detected pulsed 
7-ray emission from 11 millisecond pulsars (Abdo et 
al. 2010a, 2009a,b; Saz Parkinson et al. 2010; Guille- 
mot et al. 2011). Furthermore, the detection of radio 
millisecond pulsars (MSPs) associated with over 20 
unidentified Fermi point sources (e.g. Ransom et al. 
2011; Keith et al. 2011) has been reported, suggesting 
that the millisecond pulsar, as well as the canonical 
pulsar is one of the major Galactic 7-ray source. The 
7-ray emission from pulsars have been discussed in 
the context of a polar cap accelerator (Ruderman & 
Sutherland 1975), a slot gap (Harding & Muslimov 
2011) and an outer gap accelerator (Cheng, Ho & Ru- 
derman 1986; Takata, Wang & Cheng 2010b). 

The high quality data measured by the Fermi en- 
able us to perform a detail study for population of 
the 7-ray pulsars. For example, Takata et al. (2010b) 
and Takata, Wang and Cheng (201 Id) have studied 
the possible relation between the high-energy emis- 
sion properties and pulsar spin down power in the 
context of the outer gap accelerator model. They pro- 
posed that the outer gap accelerator model controlled 
by the magnetic pair-creation process can explain the 
observed population statistics better than that con- 
trolled by the photon-photon pair-creation process 
(e.g. Zhang & Cheng 2003). Story, Gonthier & Hard- 
ing (2007) studied the population of 7-ray millisecond 
pulsars within the context of the slot gap accelerator 
model, and predicted the Fermi observations. They 
predicted that the Fermi will detect 12 radio- loud 
and 33-40 radio-quiet 7-ray millisecond pulsars. With 
the Monte-Calro simulation of the outer gap, Takata, 
Wang & Cheng (2011a, b) have explained the observed 
distributions of the characteristics of the 7-ray pul- 
sars detected by the Fermi with the six-month long 
observation. The population studies (e.g. Kaaret & 



Philip 1996; Takata et al. 2011b, c) have also pointed 
out that unidentified MSPs located at high-Galactic 
latitudes will associate with the Fermi unidentified 
sources (Abdo et al. 2010b). 

In this proceeding, we will review our recent Monte- 
Carlo studies for the Galactic population of the 7- 
ray emitting MSPs (Takata et al. 2010b; Takata et 
al. 2011a, b,c) and the possible association with the 
Fermi unidentified sources. 



2. 7-ray emission from outer gap 
2.1. 7-ray luminosity 

For the outer gap model, the luminosity of the 7-ray 
emissions is typically 



/ Lsd, 



(1) 



where Lsd is the spin down power of the pulsar and 
the gap fraction / is defined as the ratio of the gap 
thickness at the light cylinder to the light cylinder ra- 
dius Ric = Pc/2'K, where P is pulsar rotation period. 

Zhang & Cheng (2003) have argued a self-consistent 
outer gap model controlled by the photon-photon 
pair-creation process between the curvature photons 
and the X-rays from the stellar surface. They esti- 
mated the gap fraction for the MSPs as 

fp = 7.0 X 10-2(P/i ms)26/2l(s^/lo8G)-4/^(5r5'/^ 

(2) 

where Bs is the stellar magnetic field of the global 
dipole field and 5r^ is the distance (in units of 10^ cm) 
from the stellar surface to the position, where the local 
multiple magnetic field, which dominates the global 
dipole field, is comparable to the dipole field, and it 
will be (Jrs ~ 1 — 10 cm. 

Takata et al. (2010b) argued that the incoming par- 
ticles emit photons with an energy TOgC^/o!/ ~ 70MeV 
by curvature radiation near the stellar surface and 
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that these photons can become pairs via the magnetic 
pair creation process. For a simple dipole field struc- 
ture, all pairs move inward and cannot affect the outer 
gap accelerator. However if the local field lines near 
the surface are bent sideward due to the strong mul- 
tipole field, the pairs created in these local magnetic 
field lines can have an angle greater than 90°, which 
results in an outgoing flow of pairs. In this model, the 
fractional gap thickness in this circumstance is 



0.025KP_ 



1/2 



-3 ' 



(3) 



where K ^ ^^^2*7 characterizes the local parame- 
ters. Here, -B,„.i2 and S7 are the strength of the local 
magnetic field in units of lO^^G and the local curva- 
ture radius in units of let'em, respectively, near the 
stellar surface. 

The 7-ray luminosity ([1]) can be cast in terms of 
the spin down power Lgd — 2(27r)^/i/(3c^P'*) as 



- lO''L]/^,W^6r erg s 



(4) 



for the outer gap controlled by the photon-photon 
pair-creation process, and 



3/4 _i 

8 erg s , 



(5) 
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Eq. (5); <K/Bf>~1 ■ 



Lsd(103^ergs-i) 

Figure 1: The 7-ray luminosity and the spin down power 
for the MSPs. The sole and dashed line are the relation 
predicted by the outer gap model controlled by the 
photon-photon pair-creation process and by the magnetic 
pair-creation process, respectively. The observed 
luminosity is L"*" — Aivd^F with d being the distance and 
F the observed flux. See Takata et al. 2011d. 



by the magnetic pair-creation process. Here, Lsd 34 = 

(L,,/10^^ erg s-i), K, = K/10 and Bs = BJ10» G . 2.2.1 . Charge density distribution 



In Figures [T] the model predictions given by equa- 
tions (HI) and ^ are plotted with the solid line or 
dashed line, respectively. The filled circles represent 
the MSPs detected by the Fermi-LAT. Notwithstand- 
ing the large observational errors, the data points at 
large in Figures [1] may suggest that the magnetic 
pair-creation model for the gap closing process is pre- 
ferred over the photon-photon pair-creation model for 
the — Lsd relation. 



2.2. Two layer outer gap model 

A simple description of the electric field structure 
inside the two-layer outer gap is discussed in Wang, 
Takata & Cheng (2010). We illustrate the picture 
for the gap structure in Figure [2l where we denote z 
as the coordinate perpendicular to the magnetic field 
line in the poloidal plane, and p is the charge density. 
In the main acceleration region, the charge density 
is ^10% of the Goldreich- Julian value, and an elec- 
tric field along the magnetic field line accelerates the 
electrons and positrons up to the Lorentz factor of 
r ^ 10^'^ and the accelerated particles emit ^ GeV 
photons via the curvature radiation process. In the 
screening region, the large number of pairs created by 
the pair-creation process starts to screen out the gap 
electric field. 



The previous electrodynamics studies expect that 
the particle number density increases exponentially 
near the boundary (z — hi) between the main accel- 
eration and screening regions (Cheng et al. 1986), we 
approximately describe the distribution of the charge 
density in the z-direction as follows (Fig. [2]) 



p{r) 



pi{x,(j)), if < 2 < hi{x,(j)] 
P2ix,(j}), if hi{x,4>) < z < h2{3 



(6) 



where x and (j) represent the coordinates along the 
magnetic field and the azimuthal direction. In ad- 
dition, z — and 2; = /12 correspond to the last- 
open field line and the upper boundary of the gap, 
respectively. For simplicity, we assume that the upper 
boundary (^,2) as well as the boundary (hi) is defined 
with a given magnetic field line. Because the charge 
density in the screening region will be proportional to 
the Goldreich- Julian charge density, we approximate 
that p - pGj ^ g{z, (j))pGj{r), where 



-gi{(f>), if < 2 < hi{x,4)), 
g2['^), if /ii(a;,0) < 2 < ^.2(2;, 0). 

(7) 

and gi > Q and 172 > so that \p\ < \pgj \ for the main 
acceleration region and \p\ > \pgj\ for the screening 
region. In this paper, by neglecting the 2-dependence 
of the Goldreich- Julian charge density, we approxi- 
mate as pgj{x,4>) ~ —QBx/2ttcRc, where Q and Rc 
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Figure 2: The simplified distribution of tlie charge 
density (solid line) and the corresponding accelerating 
potential (dotdashed line) of the two-layer outer gap. See 
Wang et al. 2011. 



are the angular frequency of the pulsar and the cur- 
vature radius of the field line, respectively. 

2.2.2. Accelerating electric field 

To obtain the typical strength of the electric field in 
the gap, we find the solution of the Poisson equation 
for each azimuthal angle (Wang et al. 2010), 

^$'(0;, z,0)|0=fixcd = -4:Tr[p{x, z, (j))~pGjir)]4,=Rxcd, 

(8) 

where is the electric potential of the accelerating 
field. The boundary conditions on the lower (z — 0) 
and upper (z = /12) boundaries are given by $'(a;, z = 
0, 0) = $'(a;, z = h2,4>) = 0. Imposing $' and d^' /dz 
are continuous at the boundary z ^ hi, we eventually 
obtain the solution of the accelerating electric field 
(£^11 = -d^'/dx) as 



E\\(r 



D.Bhl{x, 



-<7i(0)z'2+Ci(r>'] 



(9) 



for < z' < hi/h2 and 



cRs 



[g2mz''-l) + D,{r){z'-l)] 

(10) 



for hi/h2 < z' <\, where z' = z/h2, 

Ci{x,(f)) = -[gihi{hi - 2/12) +52(^11 - 
and 

D2{x,<j)) = -{gihj + g2hl)/hl 



In addition, we used the relations of the dipole field 
that d{Bhl)/dx - 0, dz'/dx - 0, d{hi/h2)/dx - 0, 
and approximated that dRc/dx ~ 0. 

On the upper boundary, the total potential field (co- 
rotational -|- non co-rotational potentials) in the gap is 
continuously connected to the co-rotational potential 



field outside the gap, giving d<^' /dz\z=h2 = 
h2) = 0, yielding the relation between (hi, 
(ffi, 92) as 

{h2/hif = l + (5i/.g2). 



-E 
h2 



Xz = 
I and 



(11) 



The gap thickness /12 is calculated that /12 ~ fRic, 
where / is the fractional gap thickness given by Eq. ([2]) 
or Eq. ([3]). In this paper, we present the results us- 
ing the ratio of /ii//i2 = 0.95 and the dimensionless 
charge density of main acceleration region \ ~ gi = 
0.3 (Takata et al. 2011c). The normalized charge den- 
sity of the screening region, (72 is calculated from the 
equation (fTTj) . The Lorentz factor of the accelerated 
particles is obtained by the force balance between the 
electric field and the curvature radiation drag force as 



r = 



3^ 

2e 



1/4 



(12) 



To examine the dependency of 7-ray emission prop- 
erties on the viewing geometry, we explore the typical 
3-D geometry of the outer gap (Takata et al. 2011c). 
We calculate the curvature radiation process for each 
emitting point, at which the curvature photons are 
emitted in the direction of the particle motion. Fig- 
ure [3] shows the 7-ray flux (> 100 MeV) as a func- 
tion of the viewing angle C, and of the inclination an- 
gle a. Here we calculated the rotation period from 
= 3 X 10^2 G and the gap fraction f = fp = 0.25. 
We can see that the calculated flux tends to decrease 
as the line of sight approaches to the rotation axis, 
where ^ = 0°. This is because a large viewing an- 
gle (C — J> 90°) can encounter more intense emission 
region, whereas the small viewing angle {( <^ 90°) 
will encounter the less intense region or even miss the 
emission region. Also we find in Figure [Sj a larger 
inclination angle shows a larger observed flux for a 
smaller Earth viewing angle. These dependences of 
the 7-ray flux on the viewing geometry implies that 
the Fermi has likely detected a greater number of 
the pulsars with larger inclination angles and larger 
viewing angles near 90°. 



3. Monte-Carlo simulation 

We simulate that Galactic distribution of MSPs us- 
ing the Monte-Carlo method developed by the pre- 
vious studies (e.g. Sturner & Dermer 1996). We 
would like to remark that for MSPs, the Galactic 
distribution will not depend on the spin down age 
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Viewing angie { (degrees) 



Figure 3: Dependency of the observed 7-ray flux on the 
viewing angle (^) and the incUnation angle (a). 
f = fp = 0.25. See Takata et al. (2011c). 
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Figure 4: The inclination axis a and the viewing angle ^ 
for the simulated 7-ray MSPs. Filled circles: Radio 
selected pulsars. Boxes: 7-ray-selected pulsars. See 
Takata et al. (2011c). 



of the MSPs. With the typical velocity of the ob- 
served MSPs , ~ 70 km/s, it is expected that the 
displacement of the MSPs (or binary system) with 
the typical age, > 100 Myr, becomes larger than the 
size of the Galaxy. With the relatively slow velocity, 
V < 100 km/s, however, the MSPs remain bound in 
the Galactic potential and hence their Galactic distri- 
bution does not depend on the spin down age. Because 
the independence of the Galactic distribution on the 
spin down age, we assign the "current" pulsar prop- 
erties for simulated MSPs; namely, we (1) randomly 
select the age of the simulated MSP up to 10 Gyr, (2) 
shifts the simulated MSP from its birth location to 
the current location, and (3) assign the parameters of 
the MSP following the observed distributions. 

For each simulated MSPs, we assume the radio cone 
is centered on the magnetic axis with a width de- 
scribed by the model of Kramer & Xilouris (2000). 
Because the 7-ray flux depends on the viewing geom- 
etry (Figure IS]), we assume that the viewing angle and 
the inclination angle are randomly distributed for each 
simulated MSP. 

For the sensitivity of the Fermi observations, we 
refer the sensitivity of Fermi the six-month long ob- 
servations (Abdo et al. 2010b) for the radio selected 
MSPs. However, there are no detections of the 7-ray- 
selected MSPs so far, we cannot simulate the Fermi 
sensitivity of the blind search of the MSPs. In this 
paper, therefore, we simulate the population of the 7- 
ray-selected MSPs with the Fermi sensitivity of the 
blind search of canonical pulsar. To simulate a longer 
Fermi observation, we scale the sensitivity as cx VT, 
where T is the length of the observation time. 




33 34 35 36 37 38 39 



logioLsa[erg/s] 

Figure 5: The 7-ray luminosity versus the spin down 
power for 7-ray emitting MSPs. The 200 samples of the 
simulated pulsars are plotted in the sub-figure. See 
Takata et al. (2011c). 



4. Results 

We present the results of the Monte-Carlo simula- 
tion within context of the outer gap model controlled 
by the magnetic pair-creation process, that is, the 
fractional gap thickness is given by equation [31 
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six- 


months five-years ten-years 


10° 


MSPs 


Nr 




iV, Ng Nr Ng 




Ra. Sen. (x 1) 


10 


52 


14 200 16 284 




Ra. Sen. (x 2) 


16 


48 


26 190 29 274 


10' 


Ra. Sen. (x 10) 


45 


32 


82 152 94 227 




Beaming 


106 


11 


321 41 438 62 





Table I Population of simulated radio-selected (Nr) and 
7-ray-selected (Ng) MSPs for six-month, five-year and 
ten-year Fermi observations. The first line; the results 
for ten radio surveys listed in table 1 of Takata et al. 
(2011a). The second and third lines are the results with 
the sensitivities increased by a factor of two and ten, 
respectively, and the bottom is the populations 
associated with only beaming effects of the radio 
emission. See Takata et al. 2011c. 



4.1. vs. L 



sd 



Figure m plots the distribution of the inclination an- 
gle (a) and the viewing angle (^) for the simulated 
7-ray MSPs. The fiUed-circles and the boxes repre- 
sent the radio-selected and 7-ray-selected 7-ray MSPs, 
respectively. As we can seen in Figure 21 no 7-ray 
pulsars are detected with a smaller inclination angles 
and a smaller viewing angles. This is because the 
7-ray flux decreases as the viewing angle and/or the 
inclination angle decrease (c.f. Figure [31). Hence, our 
simulation result predicts that a greater number of 
MSPs with larger inclination and larger viewing an- 
gles 90°) have been likely detected by the Fermi 
observations. 

Figuresinishows the relation between the 7-ray lumi- 
nosity and the spin down power. In the figure, we plot 
the Fermi data with errors taken from Abdo et al. 
(2010a) and Saz Parkinson et al. (2010), and present 
the simulated 200 pulsars in the sub-figures (except for 
the 7-ray-selected MSPs). The sohd and dashed his- 
tograms represent the distributions for the simulated 
and observed 7-ray pulsars, respectively. The present 
model predicts that most 7-ray canonical pulsars have 
a spin down power of Lgd ^ iQ^^-ss gj-g/g ^nd a 7- 
ray luminosity of ^ 10^^^"^^, while MSPs have a 
Lsd '-^ 10^3-35 gj-g/g and ^ j^g^^ '^^^'*-^ erg/s, which 
are consistent with the observations. 

In Figure [51 we may see that the spin down power 
Lsd and the 7-ray luminosity L-y of the simulated 
pulsars can be related as oc Lf^ with /3 ^ 
for Lsd > 1035-36 grg/g and /3 - 0.5 for Lsd < 
2^035-36 QYg/s. Because the 7-ray luminosity is pro- 
portional to L~f cx f^Lsd [Eq. ((11)], the change of the 
slope implies the switching gap closure process be- 
tween the photon-photon pair-creation process, which 
predicts oc l],^J'^ [Eq. 1^] and the magnetic pair- 
creation process, which predicts Lj oc L'^J^^ [Eq. 
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Figure 6: The distributions of Galactic latitudes of the 
Fermi unidentified sources with V < 23.21 and C > 5 
(solid line), the simulated 7-ray MSPs (dashed line) and 
canonical pulsars (dotted line). See Takata et al. (2011b). 



4.2. Population 

In Table 1, we summarize the simulated popula- 
tion of MSPs with the Fermi six-month, five-year and 
ten-year long observations. In addition, the second 
line ("Ra. Sen. (x2)") and the third line ("Ra.Sen. 
(xlO)") in the tables show the results for the radio- 
surveys increased the sensitivities by a factor of two 
and of ten, respectively, and the fourth line ( "Beam- 
ing") shows the population associated with the only 
beaming effect of the radio emission. 

With the previous radio surveys (first line in Ta- 
ble 1), the present simulation predicts 10 of the radio- 
selected 7-ray MSPs, which is consistent with the ob- 
served number 9 of the Fermi six-month long obser- 
vation. The present model predicts that 14 (or 16) 
radio-selected 7-ray MSP pulsars will be detected by 
the Fermi with five-year (or ten-year) observations. 

We see in the first line of Table 1 that the simulated 
numbers of radio-selected 7-ray MSPs increase only 
~10 sources, respectively, over even ten-year Fermi 
observations. This implies that most presently known 
radio MSPs (^ 80) might not be discovered by the 
Fermi future observations. We also see in Table 1, 
most simulated 7-ray emitting MSPs are categorized 
as the 7-ray-selected pulsars with the previous sen- 
sitivities of the radio surveys, although the Fermi 
has not confirmed the radio-quiet MSPs. We argue 
that it may be difficult to identify radio-quiet MSPs, 
because the detection of the rotation period by the 
Fermi blind search is very difficult due to, for exam- 
ple, the effect orbital motion if the MSP is in binary 
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system. We expect that the 7-ray-selected MSPs in 
the simulation correspond to the Fermi unidentified 
sources located at higher Galactic latitudes. 

We would like to remark that even we drastically in- 
crease the radio sensitivity by a factor of 10, the num- 
ber of radio-selected MSPs detected by 10 year Fermi 
observations can increase from 16 to 94, which is still 
much less than the expected 227 7-ray-selected MSPs. 
Unless the Fermi sensitivity of the blind search is im- 
proved, therefore, the most 7-ray MSPs will not be 
identified and will contribute to the Fermi unidenti- 
fied sources and/or the 7-ray background radiations. 

4.3. Fermi unidentified sources 

In Figure [SI we compare the distributions of the 
Galactic latitudes for the unidentified Fermi sources 
with the variability index V < 23.21 and curvature 
index C > 5 (solid line, Abdo et al. 2010b), the simu- 
lated 7-ray MSPs (dashed line), and canonical pulsars 
(dotted hne) with a flux Fy > 10^^^ erg/cm^s. 

We find in Figure [5] that the distribution of the 
simulated 7-ray MSPs is consistent with that of the 
Fermi unidentified sources. In particular, the MSPs 
can explain the distribution of the Fermi unidenti- 
fied sources located above the Galactic plane |6| > 10, 
which cannot be explained by the canonical 7-ray 
pulsars, which can mainly explain the unidentified 
sources located in the Galactic plane. Since the MSPs 
are in general older than the canonical pulsars, a 
higher fraction of the 7-ray MSPs, as compared with 
the canonical pulsars, is located at higher Galactic lat- 
itudes. Hence, we expect that 7-ray emitting MSPs 
are more plausible as the candidate for the origin 
of the majority of the Galactic Fermi unidentified 
steady sources located in high Galactic latitudes. 

Finally, several MSPs in the Black Window (B-W) 
systems have been discovered at the Fermi unidenti- 
fied sources (Ransom et al. 2011). Because the MSPs 
in the B-W systems are younger and have higher spin- 
down power, the 7-ray luminosity will be larger than 
that of the ordinary MSPs. Furthermore, because the 
7-ray emission from the outer gap is stronger in the di- 
rection perpendicular to the spin axis, Fermi is more 
likely to discover a greater number of MSPs with the 
Earth viewing angle ^ 90° measured from the rota- 
tion axis (Takata, Cheng & Taam 2010a; Takata et 
al. 2011c). If the angular momentum transferred from 
the accreting matter to the neutron star in the accret- 
ing stage produces the pulsar's spin axis perpendicular 
to the orbital plane, the 7-ray emissions from MSPs 
in B-W system will be stronger on the orbital plane. 
Hence, the Fermi will find the B-W systems with the 
Earth viewing angle described by edge-on rather than 
by face-on with respect to the orbital plane. In such 
a case, a greater number of the Fermi Black Window 
systems show the eclipse of the radio emissions by the 



matter ejected from the companion star. 

In summary, we have studied the population of the 
7-ray emitting MSPs. We find that the observed pos- 
sible relation between vs. L^d can be explained by 
the outer gap model controlled by the magnetic pair- 
creation process. Our Monte-Carlo study implies that 
the Fermi has detected a greater number of MSPs 
with larger inclination and larger viewing angles. Fur- 
thermore, we expect that most of 7-ray emitting MSPs 
has been missed by the Fermi observations, and still 
remain as the unidentified sources. It is likely that 
the Galactic Fermi unidentified sources located high 
galactic latitude are mainly associated with MSPs. 

We thank A.H.Kong, C.Y.Hui, P.H.T.Tam, 
R.H. H.Huang, Lupin-C.C.Lin, M.Ruderman, 
S.Shibata and R.Taam for the useful discussions. 
K.S.C. are supported by a GRF grant of Hong Kong 
Government under HKU700911P. 
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